home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / slansp.z / slansp
Encoding:
Text File  |  2002-10-03  |  4.0 KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAANNNNSSSSPPPP((((3333SSSS))))                                                          SSSSLLLLAAAANNNNSSSSPPPP((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLANSP - return the value of the one norm, or the Frobenius norm, or the
  10.      infinity norm, or the element of largest absolute value of a real
  11.      symmetric matrix A, supplied in packed form
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      REAL FUNCTION SLANSP( NORM, UPLO, N, AP, WORK )
  15.  
  16.          CHARACTER NORM, UPLO
  17.  
  18.          INTEGER   N
  19.  
  20.          REAL      AP( * ), WORK( * )
  21.  
  22. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  23.      These routines are part of the SCSL Scientific Library and can be loaded
  24.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  25.      directs the linker to use the multi-processor version of the library.
  26.  
  27.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  28.      4 bytes (32 bits). Another version of SCSL is available in which integers
  29.      are 8 bytes (64 bits).  This version allows the user access to larger
  30.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  31.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  32.      only one of the two versions; 4-byte integer and 8-byte integer library
  33.      calls cannot be mixed.
  34.  
  35. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  36.      SLANSP returns the value of the one norm, or the Frobenius norm, or the
  37.      infinity norm, or the element of largest absolute value of a real
  38.      symmetric matrix A, supplied in packed form.
  39.  
  40. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  41.      SLANSP returns the value
  42.  
  43.         SLANSP = ( max(abs(A(i,j))), NORM = 'M' or 'm'
  44.                  (
  45.                  ( norm1(A),         NORM = '1', 'O' or 'o'
  46.                  (
  47.                  ( normI(A),         NORM = 'I' or 'i'
  48.                  (
  49.                  ( normF(A),         NORM = 'F', 'f', 'E' or 'e'
  50.  
  51.      where  norm1  denotes the  one norm of a matrix (maximum column sum),
  52.      normI  denotes the  infinity norm  of a matrix  (maximum row sum) and
  53.      normF  denotes the  Frobenius norm of a matrix (square root of sum of
  54.      squares).  Note that  max(abs(A(i,j)))  is not a  matrix norm.
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSLLLLAAAANNNNSSSSPPPP((((3333SSSS))))                                                          SSSSLLLLAAAANNNNSSSSPPPP((((3333SSSS))))
  71.  
  72.  
  73.  
  74. ARGUMENTS
  75.      NORM    (input) CHARACTER*1
  76.              Specifies the value to be returned in SLANSP as described above.
  77.  
  78.      UPLO    (input) CHARACTER*1
  79.              Specifies whether the upper or lower triangular part of the
  80.              symmetric matrix A is supplied.  = 'U':  Upper triangular part of
  81.              A is supplied
  82.              = 'L':  Lower triangular part of A is supplied
  83.  
  84.      N       (input) INTEGER
  85.              The order of the matrix A.  N >= 0.  When N = 0, SLANSP is set to
  86.              zero.
  87.  
  88.      AP      (input) REAL array, dimension (N*(N+1)/2)
  89.              The upper or lower triangle of the symmetric matrix A, packed
  90.              columnwise in a linear array.  The j-th column of A is stored in
  91.              the array AP as follows:  if UPLO = 'U', AP(i + (j-1)*j/2) =
  92.              A(i,j) for 1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) =
  93.              A(i,j) for j<=i<=n.
  94.  
  95.      WORK    (workspace) REAL array, dimension (LWORK),
  96.              where LWORK >= N when NORM = 'I' or '1' or 'O'; otherwise, WORK
  97.              is not referenced.
  98.  
  99. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  100.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  101.  
  102.      This man page is available only online.
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.